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Abstract 



Three-wave interactions (M <-> P + Q) of plasma waves propagating parallel to the 
mean magnetic field at frequencies below the electron cyclotron frequency are con- 
sidered. We consider Alfven-ion-cyclotron waves (L), fast-magnetosonic-whistler 
waves (R), and ion-sound waves (I). Especially the weakly turbulent low-beta plas- 
mas like the solar corona are studied, using the cold-plasma dispersion relation for 
the transverse waves (L, R) and the fluid-description of the warm plasma for the 
longitudinal waves (I). We analyse the resonance conditions for the wave frequencies 
uj and wavenumbers k (i.e., wm = oop + uq and ku = kj> ± &q, and the interaction 
rates u MPCi of the waves for all possible combinations of the three wave modes, and 
list those reactions that are not forbidden. One of the waves has to be longitu- 
dinal and two transverse. This demonstrates an attractive feature of the theory: 
the conservation of angular momentum is implicitly built in. In a low-beta plasma, 
non-zero reaction rates are obtained for (i) L + <-> I + + L~ and R + <-> I + + R~ in a 
wide frequency range extending from the MHD frequency range to the resonances 
of the waves; (ii) for L + <-> L + + I + in more narrow frequency range, where at least 
one of the L waves is in the dispersive frequency range; and (iii) I + <-> L + + 
close to the resonance of the L mode. The reaction types (ii) and (iii) have, to our 
knowledge, not been discussed before in low-beta plasmas. In high-beta plasmas, 
reactions I + <-> R + + L~ and I + <-> L + + R are the main reactions, extending 
down to the MHD frequencies and discussed earlier for the non-dispersive case, but 
new reactions involving dispersive waves are found in high-beta plasmas as well: 
I + <-> L + + R + and R + <-> I + + ^ are now possible in limited frequency ranges 
involving at least one dispersive transverse wave. We discuss the implications of the 
discovered new reactions to turbulent cascading in space plasmas. 
Keywords: plasma waves - turbulence - wave-wave interactions 



1 Introduction 



Three-wave interactions constitute the lowest-order non-linear coupling between 
wave-modes in plasmas. These interactions are of form, where either two modes 
coalesce into one or one mode decays into two, i.e., M <-> P + Q, with momentum and 
energy being conserved at microscopic level, i.e., wm = uj p + ujq and k M = k P + 
Thus, three-wave interactions provide the leading-order description of the evolution 
of weak turbulence, where it is assumed that the amplitudes of the interacting wave 
packets are small enough for a perturbation theory approach to apply. Turbulence 
evolution, on the other hand, is one of the key issues in plasma astrophysics and 
space physics. It is intimately related to fundamental questions like plasma heating 
and particle acceleration and transport in collisionless plasmas. To get started with 
weak turbulence theory, the reader is referred to. 1 

In this paper we will investigate three-wave interactions of plasma waves prop- 
agating parallel to the mean magnetic field in a magnetised, low-beta plasma. The 
pioneering work by Chin and Wentzel 2 considered three-wave interactions between 
low-frequency plasma waves under the random-phase approximation, assuming that 
the waves fulfill the MHD dispersion relations, i.e., uj = vp^k\\ for the transverse 
waves and uj = vsk\\ for the longitudinal waves. Their study revealed that for 
a low-beta plasma, three-wave interactions where an Alfven wave decays into a 
parallel-propagating sound wave and an anti-parallel-propagating Alfven wave, i.e., 
A + <-> S + + A~, dominate the evolution of weak turbulence. In these interactions, 
the number of Alfven wave quanta is conserved, and since the sound waves are 
generally easily damped, the interactions proceed typically in the direction of wave 
decay. Thus, their study predicted an inverse cascade of energy from high to low fre- 
quencies with a simultaneous heating of the plasma by the dissipation of the sound 
waves emitted by the three-wave interaction. Their model has been applied by nu- 
merous authors to study different astrophysical scenarios from the heating of the 
solar corona 3 to the confinement and transport of cosmic rays (CR) in the galaxy 4 
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and more recently to CR acceleration in astrophysical shocks. 5 The work of Chin 
and Wentzel, however, was developed under the MHD approximation. In this work, 
we will consider an extension to the work by adopting dispersion relations more 
suitable for a collisionless plasma. 1,6 Specifically, we will consider interactions be- 
tween left- and right-hand circularly polarised transverse waves and longitudinal ion 
sound waves at frequencies well below the electron gyro-frequency. We will apply 
the results of the analysis to conditions in the solar corona, pointing out several new 
evolution scenarios of parallel propagating plasma waves, which have relevance to, 
e.g., coronal heating and particle acceleration at the Sun. This treatment is con- 
sidered a kinetic treatment for two major reasons: First of all we are not using the 
limited MHD dispersion relations which would not hold for frequencies well above 
the ion gyrofrequency. Secondly we are not bound to thermal distribution functions 
of the background plasma, which means we do not depend on collisions to produce 
thermal distributions. As will be shown later in this paper for Maxwellian distri- 
butions, however, the actual interactions only depend on the parameter (3 and the 
ratio Qi/ujpi. 

The structure of the paper is as follows: in §2 we will summarise the theory of 
three-wave interactions for collisionless plasmas, in §3 we will analyse in detail the 
three-wave interactions between plasma waves propagating parallel to the magnetic 
field, in §4 we will discuss the results and present the conclusions from the analysis. 

2 Three- wave interactions in a collisionless plasma 

The framework used here to describe the interaction of waves in collisionless plasmas 
employs a quasi-particle description of waves derived previously by 2 and described 
in detail by 1 and references therein. These theories are, as already mentioned above, 
usually derived for a fluid plasma (for 2 ) or an unmagnetised plasma (as in 1 ). 

Our calculations in the kinetic limit for magnetised plasmas follow 7 using the 
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interaction tensor from. 8 The physical picture is the following: The interaction 
of waves is transmitted by the electric field of these waves. From this interaction 
an infinite series of interaction terms is given, where the linear (first order) term 
describes the propagation and the first nonlinear (quadratic) term is the simplest 
form of interaction. 

Given three reservoirs of wave quanta denoted by N M , N p , and the evolution 
equation for the occupation number of the wave quanta is given by 

f d3fc/ / d? k " 7 .mpq/» y u»\ 
J (2nfJ (2vr)3 M 

x {N p (k')N Q (k") - N M (k)[N p (k') + 7V Q (fc")]} (1) 

/ d3/c f d3/c " ..mpQ/l. y u»\ 
J (2tt)3 J (2vr)3 M ^ ' 

x {N p (k')N Q (k") - N M (k)[N p (k') + iV Q (fc")]} (2) 

r _d^k_ r d 3 k' mpq/1 ./ 
./ (2tt)3 y (2vr)3 M ^ ^ 

x {7V p (fc')^ Q (fc") - 7V M (fc)[^ P (fc') + 7V Q (fc")]} ■ (3) 

Here, the primed [double-primed] quantities refer to the mode P[Q] and the unprimed 
quantities to mode M. The interaction rate is defined as 

8(2n) 7 he 4 W| Wj_ 
u) M u p u<* W p W% X 

\K MPQ (k, k', k")\ 2 5 3 (k - k! - k")8{u -J - lu") (4) 

The interaction can be derived when using the vector potential of the electric field 
(in temporal gauge). From the vector potential of all involved waves the response 
tensor can be calculated, which is basically given by the quantity u MPCl . In weak 
turbulence theory the response tensor is calculated only to the second order in the 
vector potentials. The fact that the electric fields are interacting is resembled by 
the factor We/Wt, which simply describes the ratio of electric to total energy. The 



dN M (k) 
dt 

dN p (k') 
dt 

dN®(k") 
dt 



u MPQ (k,k',k") = 
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physical details of the response tensor are now hidden in n MPCi : 



K MPQ (fc, fc', fc") = (eJ 4 )*eJe?K ij -,(fe, fe', w'; fe", w") 
1 



K>ijl(k : fc , fc ) 



fe , fe ) + Knj(k, k , fc ) 



(5) 
(6) 



The quantity /t MPQ is, thus, a scalar formed by contracting the three polarisation 
vectors of the interacting waves with the quadratic response tensor of the plasma. 
For our further calculations we will make use of the quantities k and e, only. The 
form of the k tensor is given by 8 



K 



iji 



<Wv+^'^ir 27r / d P\\ I P±dp±e 



+oo 



3 ie(pV'-/nV) 



X 



X 



X 



LU — fJ,Q — k\\V\\ 

a* (J, fc', p) A ^ ~ V ' ] , fc', p) + 7 *(u/, fc', p) A 
J <9p ± p ± J J dp\\ 



X 



cj" -iAI- k"v\\ 



a*{u",k",v)^-+^{u\k",v)^- 
op± op\\ 



fip±,p\\) 



(7) 



Here, / is the distribution function in momentum space, n is the particle density 
of the background medium, and q is the particle charge. We will limit ourselves to the 
case of gyrotropic non-relativistic /, since relativistic or non-gyrotropic distribution 
functions will increase the complexity of the problem drastically. In the derivation of 
the equations for the interaction rates no further assumptions are made. The general 
form of the interaction tensor does not depend on the specific form of the distribution 
function. It is, however, necessary to choose a specific form of the distribution 
function to determine the possible interactions (which we do in this paper) and to 
compute the actual interaction rates from the derived general equations (which we 
leave for the future). We use the dispersion relation of ion sound waves (and partially 
also those of L- and R-waves) derived from a Maxwellian distribution function. Using 
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another form of the distribution could, in principle, change the possible interactions 
if the resulting dispersion relations are qualitatively different from those obtained 
from the Maxwellian one (e.g., if the assumed monotonic behavior of u(k) would be 
violated). 

The summation of //, p, u, u' , v" is going from — oo to oo. ip is the phase angle of 
the wave vector around the magnetic field, the velocity V is defined as 



^[e^J s -i(2)+e-^J s+1 (z)] ' 



V(s,k,v) 



v\\J s (z) 



J 



(8) 



and the remaining coefficients as 



oti(u, k, s 
7i(w, k,s 
k, s 

P 2 {uJ,k,s 

P 3 (uj,k,s 



Vi(s, k,v) - 5 i3 (uj - sQ - k\\v\\) — J s (z) 

v± v± 

k\\Vi(s, k, v) + 5 i3 (uo — sVL — k\ i \v\\)J s (z) 
v± 

IQ — k\\V\\ 



v± 
v±k± 

v± 

\q\ 

k±v± 



-V 2 (s, k, v) + k±v± smi/j J s (z) 
Vi(s, k, v) — k±v± cosip J s {z) 



smip Vi(s, k, v) — cosip V 2 (s, k, v)] 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 

(15) 



When we restrict the discussion to parallel propagating waves, i.e., k± = 0, we obtain 
J s (z) — > 5 s0 and all terms involving the phase angles of the waves vanish. Thus, we 
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get 



i 



where 



E 



X 



c J-oo Jo v-nll-k\\v\i 



X 



&) (a/, Mi , p) ^ - yf 3 Ml > + 7? Ml > p) 4 

^(^Mi'^^ + ^^.M^)^ 



x 



(16) 



V(s,kn,v) 



7*(w,fc||,s) 



7 ^(5 S)1 + 5 S ,_1) 



\ V ll^,0 



/ 



a* (a;, fcii,s) = 



(3 (u,kn,s) = 



k\\V\\ 



5 a ,i + 5 s ,-i 
ie(6 8>1 - 6a-!) 

V 

^ie(8 St i - 5 S _i)^ 


\ 



V 



1 + 

ie^(5 S) i - 5 s _i) 



(17) 



(18) 



(19) 



(20) 



We will limit ourselves to propagation of waves parallel and anti-parallel to the 
magnetic field. In this case the resonance conditions for any three-wave interaction 
M <-> P + Q read 



k = k'±k" 



uj = J + cu" 



(21) 
(22) 



where the wave numbers, k = \k\ etc., i.e., the lengths of the wave vectors || £3, and 
the frequencies are taken to be positive. Here, as above, the primed [double-primed] 
quantities refer to the mode P[Q] and the unprimed quantities to mode M. 



3 Results 

3.1 Interaction rates 

The wave modes propagating parallel to the magnetic field can be either purely 
transverse (T) or purely longitudinal (Lo). In the first step we will now show, which 
of the interactions between these wave modes are forbidden because of vanishing 
interaction rates. This can be done by finding the vanishing components of the k 
tensor and the symmetries of the k tensor, which cause the vanishing contraction of 
the polarisation tensor E = e M e p e^ (components = efeje^ ) with the n tensor. 
The results for all possible three-wave reactions of parallel propagating waves, to be 
discussed in detail below, are summarised in Table [Tj 

The components of the R, tensor itself vanish whenever none or two of the indices 
are equal to 3 



(23) 

(24) 

(25) 

(26) 



Kill 


= K112 


= «121 


= ^122 


«211 


= ^212 


= ^221 


= ^222 




^133 


= «313 


= ^331 




^233 


= ^323 


= ^332 



When we take into account that the polarisation vectors of transverse waves only 
have non-zero 1 and 2 components, while the longitudinal waves have a non-zero 
3 component, we can make the first statement concerning the wave modes that 
may interact. Assuming interactions between three transverse wave modes, the 
polarisation tensor's components are zero, whenever one index is 3. Multiplying 
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Table 1: Three-wave interactions between waves with wave vectors parallel to the 
ambient magnetic field. Specific dispersion relations analysed only at frequencies 
below the electron cyclotron frequency. 





\/ Q 1 1 rl in 
V dllU. Ill 


Ullllllcll L o 




QjJlll UUllot:! VcU. : 


T^T + T 




= for all non-zero E^ 




No 


Lo <-> T + Lo 


— 


Kiji = for all non-zero Eyi 




No 


L <-> R + I 


— 


terms cancel to produce EijiKiji 


= 


No 


R^ L + I 




terms cancel to produce E^k^i 


= 


No 


I <-> L + L 




terms cancel to produce E^k^i 


= 


No 


I <-> R + R 




terms cancel to produce E^k^i 


= 


No 


L+ <-> L+ + I- 




requires dui^/dlc < 




Yes 


R + <-> R+ + I- 




requires du^/clk < 




Yes 


R+ <-> 1+ + R" 


low j3 


allowed in MHD 




Yes 




high /3 








L+ <-> 1+ + L" 


low /3 


allowed in MHD 




Yes 


1+ <-> L+ + R- 


high (3 


allowed in MHD 




Yes 




low /3 








1+ R + + L" 


high (3 


allowed in MHD 




Yes 


1+ <-> L+ + R+ 


low /3 






Yes 




high /3 






Yes 


L+ <-> L+ + 1+ 


low (5 


decaying cj l < 




Yes 


R + <-> R + + 1+ 


high /3 


decaying cj r ~ f2j 




Yes 
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this with the k, tensor results in a zero. So interactions between three transverse 
wave modes will not take place. The same applies for interactions involving two 
longitudinal and one transverse wave. 

To make further assertions about possible interactions we shall now investigate 
the symmetry of the k, tensor. The index symmetries, 



«113 = K 223 (27) 

K 123 — — K 213 (28) 

«131 = K 232 (29) 

«132 = -«231 (30) 

«311 = K 322 (31) 

«312 = -«321, (32) 



can be calculated either directly or simply by noting that the tensor is necessarily 
transversely isotropic, i.e., that there has to be a rotational symmetry around the 
rr 3 -axis aligned with the magnetic field and all the wavevectors, as this is the only 
preferred direction in the system. 

Before we can determine the components of the polarisation tensor, we have to 
specify which wave modes are used. We will now limit our discussion to the L- and 
R-mode for the transverse wave fulfilling the dispersion relations 

N 2 — L — 1 — . ^— - . ^g- (33) 

and 

i .2 , ,2 

N 2 = R=l- . p ' . - . pc , 34 
^(cj + iZi) u(u — \l e ) 

respectively. Here, iV = ck/u is the refractive index of the wave and L and R 
are the Stix-parameters. 9 fi; and Q c are the gyrofrequencies of ions and electrons, 
respectively. The plasma frequencies are denoted by aj pi and u pc . The polarisation 
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Table 2: Non-zero elements of the polarisation tensor E^i = efeje z for interactions 
of the type T <->• T + Lo. 







•^113 


E\23 


-^213 


E223 


L 


«-> L + I 


1/2 


-i/2 


i/2 


1/2 


L 


<-> R + I 


-1/2 


-i/2 


-i/2 


1/2 


R 


<-> R + I 


1/2 


i/2 


-i/2 


1/2 


R 


<-> L + I 


-1/2 


i/2 


i/2 


1/2 



vectors associated with these wave modes are 



i/V2 




i/V2 

V / 



(35) 



For the longitudinal mode we choose the ion-sound mode 



U! 



2 - u 2 \l) 



(36) 



where c s = ^k-QT e /mi is the ion sound speed, k& is Boltzmann's constant, T e the 
electron temperature, m ; the ion mass and Ad the Debye length, with the polarisation 
vector 

(37) 

vv 

The non-zero components of E^i for the interaction pairs of type T <-» T + Lo 
are given in Table |2j Combining this with the k, tensor, we find that the interaction 
rate vanishes identically for the interactions L <-> R + 1 and R <-> L + 1, while for the 
other two pairs it is 

(38) 



^ocKia-i^l 2 



which is non zero. So for this type of interactions only L <-> L + 1 and R «-> R + 1 
may have non-zero interaction rates. 
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Table 3: Non-zero polarisation tensor elements for interactions of the type I <-» T+T. 







-^311 


-E312 


E321 


E322 


I «-> 


L + L 


-1/2 


i/2 


i/2 


1/2 


I «-> 


R + R 


-1/2 


-i/2 


-i/2 


1/2 


I «-> 


R + L 


-1/2 


-i/2 


i/2 


-1/2 


I «-> 


L + R 


-1/2 


-i/2 


i/2 


-1/2 



For the next set of interactions, i.e., of type I <-> T+T, we have given the non-zero 
components of the polarisation tensor in Table |3j As the symmetry of the K-tensor 
does not change, we find that here the set of allowed interactions is I ^ R + L and 
I «-> L + R. 

Reviewing the set of eliminated interactions one finds that these selection rules, 
which were derived solely from the interaction rate without explicit regard to the 
resonance conditions, are going back to angular momentum conservation. When 
one associates a spin of +H with an L-wave and —h with an R-wave, 1 it is obvious, 
that the discarded reactions are the ones with non-conserved angular momentum. 
All of the remaining interactions are, thus, conserving angular momentum at the 
microscopic level. With this reduction of possible interactions we will now examine, 
which interactions fulfill the resonance conditions. 

3.2 Resonance conditions 



Because of the resonance conditions, Eqs. (21) and (22), the interactions 



L+ «-> L+ + T (39) 
R+ <-> R+ + r, (40) 



can be discarded immediately. The transverse wave has the same dispersion relation 
on both sides, and it has duj/dk > 0. However, the frequency of the transverse wave 

1 Note, that the propagation direction of the wave does not play a role, since the plasma physics 
definition of the handedness of circular polarisation is with respect to the magnetic field. 
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should be decreasing from left to right in the reaction while its wavenumber should 
be increasing, which is not possible. Thus, although these reactions do conserve the 
angular momentum, they cannot occur in plasmas. 

The resonance conditions are next carefully evaluated for the remaining reactions. 
It is practical to perform this separately in two distinct parameter regimes: the high 
and low plasma (3 regions. 2 This distinction is useful as the phase speed of the sound 
wave is below that of the low-frequency transverse waves in the low-beta case and 
above it in the high-beta case at low frequencies. 

3.2.1 Low /3-plasma 

The dispersion relations of parallel propagating transverse and longitudinal waves 
in a low-beta plasma are depicted in Fig. [T} In the low-/? case the frequency of 
R-waves well below the electron cyclotron resonance is always higher than that of 
the longitudinal wave with the same k. This applies also for the L-wave unless it 
approaches the ion gyroresonance. 

The first two interactions to be inspected are 

L + <-> I+ + L" (41) 
R+ <-> I+ + R- (42) 

which have been discussed under the MHD approximation already in. 5 They fulfill 
the resonance condition also in the dispersive limit up to the wave number, where 
a>r(fc) ~ u^ (k), so they are allowed in the low /3-regime. We will analyse these 

reactions numerically below. 

2 In this paper we define plasma (3 to be the ratio of the squared ion-sound speed to the squared 
Alfven speed, (3 — c 2 /V^, which is about one half of the ratio of gas pressure to magnetic field 
pressure, if T c ~^> Ti. 
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Figure 1: Dispersion relation of parallel-propagating plasma waves in low-beta 
plasma, i.e., for c s /va = 0.1. The red (solid) curves represent the ion sound wave 
for a high-density (co> pi = 10 fl e ) and low-density (cj pi = 100 Q{) plasmas, the blue 
(dot-dashed) curve is the L wave and the green (dashed) curve the R wave. Note 
that the horizontal parts of the dispersion relation represent very strongly damped 
waves. 
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The interactions 



L + <-> L+ + I+ (43) 
R + «-> R+ + I+ (44) 

can not fulfill the resonance conditions in the MHD regime, since in the low-frequency 
range the phase speed of the transverse waves is larger than that of the longitudinal 
wave for low (3. The interaction involving the R-wave may work only if the decaying 
wave frequency is very close to wr = Q e , while the L-wave interaction is limited to 
work in a small frequency range close to fij. We will analyse the reactions numerically 
below. 

The L-wave reaction can be analysed analytically using the approximate disper- 
sion relations, 



W - f J^- (45) 
h(uj) » -. (46) 



The resonance conditions yield 



L) / Oj L)i / fli UJ — iO\ 



Va V — ^ Va V — ^1 



(47) 



where is the frequency of the decaying L wave and iO\ the frequency of the daughter 
L wave, uj — uj\ being the frequency of the I wave. This condition can be written in 
the form /(w/fii) = /(wi/fii), where (Fig. [§ 

/(x) = -^--^=, /3 = c s 2 /^ A 2 (48) 

VP V 1 - £ 

The function has two zeros at x = and x = 1—/?, and a positive value in betweenand 
a single maximum. The values of uj\ and uj, thus, correspond to the two roots of 
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f[x] 



0.1 



0.01 



0.001 



X 



10 



-4 



0.001 



0.01 



0.1 



1 



Figure 2: The function of f(x) used to analyse the resonance conditions in the 
interaction L + <-> L + + I + . The value of (5 used in this plot is ~ 0.0087. 

the equation f{u/Vt\) = C, where C is a constant between and the maximum of / 
which for a very low beta is at x « 1 — (^(i) 1 ^. For higher values of (3, the zero of 
the derivative, x max , can be found by iteration taking Xmlx = 1 — (t/3) 1 / 3 and 



x 



(n+l) = l 



'1„W 
, 2 x max 



1) 2 ] 1/3 , 



(49) 



Thus, one of the roots (the daughter wave oj\) is at < x < x max and the other 
one (the decaying wave u) at x max < x < 1 — Cg/V^. The difference of the roots 
(multiplied by gives This problem may be solved analytically, but the formula 
is of little practical use. 

The R-wave reaction can be analysed in a similar manner, but the details are 
left to the interested reader. 
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By similar arguments, the three interactions 



1+ <-> L+ + R+ (50) 
1+ <-> L+ + R- (51) 
1+ <-> R+ + LT (52) 

have no solution for the resonance condition in the low-frequency limit as the phase 
speed of the transverse waves exceeds the phase speed of the longitudinal wave. As 
above, however, this is no longer true in the direct vicinity of the resonances of the 
transverse waves. 

For simplicity, we will consider a high-density plasma, for which a linear dis- 
persion relation for the longitudinal wave is valid at all frequencies below fl e . The 
resonance conditions for the first two reactions yield 



UL + ^ R = k L (uj L ) ± k R (uj R ), (53) 



i.e. 



*lW-- = -T*bW. (54) 

C s C s 

The left-hand side is positive only for k\, > fli/c s . The right-hand side is always 



positive for the lower sign, interaction (51), and for the upper sign, interaction (50), 



as long as k R < f2 e /c s . Similarly, for the third interaction (52) one can conclude 
that solutions are possible only for k R > Q c /c s . Thus, the following reactions are 
possible if u pi > Q c : 

I + ^ L + + R + for k L > Vtjc s and k R < £l e /c s (55) 

I + ^ L + + R + for k L < Vljc s and k R > fi e /c s (56) 

1+ <-> L+ + R- for k L > a/c s (57) 

I + <-> R+ + L~ for k R > tt c /c s . (58) 
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Note that if uj v \ < Q e , interactions involving electron-cyclotron waves, Eqs. (56) and 



(58), are not possible at all, and the ranges of validity of the other reactions become 



modified (i.e., u R < cj pi — Qi needs to be satisfied). 



In the reactions involving the ion-cyclotron waves, Eqs. (55) and (57), we may 
write lj l = and get 

k L = ki(uj R + ni) =f k R (u R ), (59) 

which gives an analytical solution of the wave number of the left-handed wave as a 
function of the frequency of the right-handed wave. Note that this equation is valid 
also for Qi < u pi < Q c . 

In summary, from the seven distinct possible interactions not eliminated by an- 
gular momentum conservation and du/dk > 0, only two can take place over the full 
frequency range in a low /3-plasma with parallel wave propagation: 



L+ 



I+ + L- 
1+ + rV 



In addition, the interactions (43)-(44) and (50)-(52) may have solutions where one 



of the waves are in the electron- or ion-cyclotron range. Reactions (43), (50) and 



(51) all yield solutions where an ion-cyclotron wave interacts with an MHD wave. 



Numerical analysis. For the two valid interactions and the one interaction with 
limited range we will now show the frequency triads for the resonance condition. As 
a model system for a low (3 plasma we have chosen typical parameters for a coronal 
hole 



Number density ri\ = n e = 2 x 10 cm 3 
Electron temperature T e = 1 x 10 6 K 
Magnetic field B = 2 G 
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The values correspond to (3 = c 2 s /Vl « 0.0087 and u; pi /fii « 310. 

In the upper right panel of Fig. [3] the interaction I + + L~ «-> L + is shown. 
This interaction is very similar to the one described in, 5 but for frequencies close to 
the ion gyrofrequency f2 ; the effect of the L-wave resonance can be seen. At those 
frequencies, therefore, an ion-cyclotron wave is able to decay into a low-frequency 
Alfven wave and a high-frequency ion sound wave. This process, however, has to 
compete against the dissipation of the mother wave due to ion cyclotron resonance. 

The interaction I + + R~ «-> R + (upper left panel of Fig. [3]) is similar to the 
previous one in the low-frequency region, but shows a different behaviour for uo > Q u 
since the R-wave has no resonance there. The R-waves (the decaying one and the 
daughter wave) have a linear frequency relation, while the ion sound wave resembles 
the non-linear frequency difference due to the non-linear dispersion relation of the 
R-wave. It is noteworthy that the inverse cascading of the R-waves may proceed all 
they way from close to the electron-cycotron frequency through the Whistler range 
to the MHD regime. 

In Fig. [4] the limited ranges of validity for the reactions L + + I + <-» L + and 
I + *-> L + + R are plotted. The dispersive features of these new interactions are 
clearly visible. 

3.2.2 High-/? Plasma 

As we have done for the low (3 case, we will also analyse the interactions with non- 
vanishing interaction coefficients for a high-beta plasma. We will, however, keep the 
discussion based on analytics more limited in this case than for the low-beta case. 
The dispersion relations in high-beta plasma are depicted in Fig. [5j 
For /3 > 1, the interaction 

1+ <-> R+ + L+ (60) 
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Low yS-Plasma 



r + r < — > r + r + l < — > l + 




High y8-Plasma 




Figure 3: The frequencies of the daughter waves as a function of the frequency of 
the decaying wave in interactions, which fulfill the resonance condition in the whole 
frequency range including the MHD regime, analysed earlier in. 2,5 In the plots L- 
waves (blue [dot-dashed] curve), R- waves (green [dashed] curve) and I- waves (red 
[solid] curve) can be seen. 
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Low /3-Plasma 
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Figure 4: The frequencies of the daughter waves as a function of the frequency of 
the decaying wave in the new interactions, which fulfill the resonance conditions 
only in the dispersive wave regime in a low-beta plasma. In the plots L-waves (blue 
[dot-dashed] curve), R- waves (green [dashed] curve) and I- waves (red [solid] curve) 
can be seen. 
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(j/fLi 

w pi = 10 fi e 




Figure 5: Dispersion relation of parallel-propagating plasma waves in high-beta 
plasma, i.e., for c s /va = 10. The red (solid) curves represent the ion sound wave 
for a high-density (u P i = 10 Q c ) and low-density (co> pi = 100 fli) plasmas, the blue 
(dot-dashed) curve is the L wave and the green (dashed) curve the R wave. 

cannot fulfill the resonance condition at low frequencies, as the phase speed of lon- 
gitudinal wave is higher than that of both the transverse waves. However, if the ion 
sound speed is below the maximum value of the whistler phase speed (~ 20 Va) or 
if the plasma has a low-enough density so that uj pi < Q c , the R-wave phase speed 
exceeds the sound-wave phase speed (Fig. [5]) and reactions may occur. 
The interaction 

L+^I+ + L- (61) 

is not possible in a high-/5 plasma, as ki,(u) — uj/c s > V uj < O; in a high- 
beta plasma, and this makes it impossible to fulfill the resonance condition. The 
interaction 



L + «-> I+ + L+ (62) 
will not work either, because ki,(u) — u/c s is a monotonically increasing function, 
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which allows the resonance condition, k^u) — uj/c s = k^(ui) — uj\/c s , to be fulfilled 
only at uj = uj — u)\ = 0. 
However, the interactions 

R+ «-> 1+ + R* (63) 

are possible as long as at least one of the R waves falls in to the dispersive frequency 
range. The resonance conditions for a linear sound wave, U\ = c 8 ki, yield 



k R (u) -- = HbrM - -. (64) 



The upper sign corresponds to g(y) = g(yi) with y = uj/VL e and (Fig. |6 



n V / m e mi 

V [m e m i + y)[l-y) 



which has zeros at y = and 



1 - m e /mi ± ^/(l - m e /mi) 2 - 4(/3 - l)m e /m- l 

y± = g ^ ) 

y_ ^ (P-\)mJm h y + ^l-/3m e /mi (67) 



and a single minimum between and y_ and a single maximum between y_ and 
y + . Thus, R waves with frequencies y above the position y max of the maximum are 
decaying into R waves with frequencies below y m3X , analogously to the low-beta case 
L + <-> L + + I + , but now the mother- wave frequency is very close to Q c . However, 
there is also another branch of interactions, which corresponds to the negative values 
of g. There, R waves with frequencies above the position y min of the minimum of g, 
i.e., with y e (y m m,y-) are decaying into R waves at frequencies y G (0,y min ). The 
value of y m \ n can be found iteratively in a similar manner as x max for f(x) above, 
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Figure 6: The function of g(y) used to analyse the resonance conditions in the 
interaction R + R + + I + . The dashed curve gives \g\ for g(y) < and the solid 
curve the positive values of g. The value of j3 used in this plot is ~ 14. 

since at y < y- <C 1 and with x = m\yjm e we have 



g 

dg 
dy 



m e x 

dxdg 1 \x + 1 



dy dx y/p (1 + xfl 2 



(68) 
(69) 



giving, for g'{y m \ n ) = 0, an iterative formula: x^J = 1 and 



X min ^ ~~ [@(2 X niL + l) 2 ]^ 1- 



(70) 



Thus, these solutions are found for mother- wave frequencies ui G [x min f2i, (f3 — 



The lower sign in Eq. (64) yields solutions 



m e /mi 



m e /mi 



fj3 \ (m e /mi + y)(l - y) </j3 V (m e /mi + yi)(l - yx) 



(71) 
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meaning that the solutions are given by the roots of 



, 1 / m e /mi . 

g(y) = jfr _ + J ^ 1 , 72 

rj3 V (m e /mi + - 2/1) 



where the right-hand side is monotonic and can have any value between and 00. 
Thus, this reaction involves mother- wave frequencies in the whole range where g(y) 
is positive. 

Finally, we are left with two interactions, which have been treated in 5 for the 
non-dispersive case 

1+ <-> L+ + R- (73) 
1+ <-> R+ + LT. (74) 

These should work also in the dispersive case. We will analyse these reactions in 
more detail numerically. 

Numerical analysis. As for the low /3-case we plot the solutions for the resonance 
condition in the high /3-case. As a model system for a high (3 plasma we have chosen 
parameters typical for the downstream region of a strong coronal shock 

Number density rij = n e = 8 x 10 9 cm -3 (75) 
Electron temperature T e = 2.5 x 10 7 K (76) 
Magnetic field B = 5 G (77) 

The value of beta for this plasma is (3 — c^/Vl ss 14 and u P i ~ 1.3fi e - 

The first interaction in the high-/5 case has also been discussed in: 5 I + <-» L + + R~ 
(lower left panel of Fig. [3]). The solution of the resonance condition in the dispersive 
case is similar to the non-dispersive case for uj < Q^. In the high-frequency regime 
we find, that the L-wave goes into resonance. For low frequencies the interaction 



24 



takes place for an L-wave of slightly higher frequency than the R-wave, this changes 
at w = H;. For high frequencies the energy is supplied almost completely by the 
R-wave, while the momentum comes from the L-wave. 

The next interaction has also been discussed in: 5 I + <-> R + + L~ (lower right 
panel of Fig. [3J. Here the dispersive effects show a completely different behaviour. 
The L-wave frequency drops to zero for high frequencies, while the R-wave shows 
an almost linear increase. 

For the interaction R + <-> R + + 1 + (top panel of Fig. [7]) it can be clearly seen that 
there is only a small frequency band where the resonance condition can be fulfilled. 
This is due to phase speed of the R-wave which exceeds that of the sound wave only 
for a limited frequency range. 

As in the previous case, the interaction R + «-> I + + R (middle panel of Fig. [7]) 
shows clear dispersive effects. It is limited to the Whistler frequency range but has 
a larger range of validity than the previous case. In contrast to the previous case 
the upper limit is now given by the electron cyclotron limit of the R-wave. 

Finally, the reaction I + <-> R + + L + is plotted in the bottom panel of Fig. [7j This 
shows that ion-cyclotron waves are able to interact with higher-frequency Whistlers 
propagating in the same direction, producing longitudinal waves that are rapidly 
damped. This reaction, thus, provides a dissipation mechanism for Whistlers in 
plasmas with cold ions and hot electrons. 

4 Discussion and conclusions 

We have analysed three-wave interactions of parallel-propagating plasma waves, con- 
centrating on the effects of wave dispersion on the interactions. Our analysis shows 
that the theory is consistent with the conservation of angular momentum at micro- 
scopic level: the total amount of spin carried by the wave quanta must be conserved 
in the three-wave interactions. The reactions conserving the angular momentum are 
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Figure 7: The frequencies of the daughter waves as a function of the frequency of 
the decaying wave in the new interactions, which fulfill the resonance conditions 
only in the dispersive wave regime in a high-beta plasma. In the plots L-waves (blue 
[dot-dashed] curve), R- waves (green [dashed] curve) and I- waves (red [solid] curve) 
can be seen. 
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further analysed by employing the dispersion relations of left- and right-hand circu- 
larly polarised waves at frequencies below Q { and Q c , respectively, and of ion-sound 
waves at frequencies below u P i. 

Our analysis extends the previously analysed 2,5 low-beta interactions, 

L + <-> I+ + L - 
R + <-> I+ + R- 

to higher frequencies and shows that the latter reaction may occur also in high-beta 
plasma at Whistler frequencies of the mother R wave, where its phase speed exceeds 
that of the ion sound wave. The implied efficient inverse cascading of the R waves 
has important consequences on energetic particle transport and acceleration in col- 
lisionless plasmas: as higher-energy particles generally interact with lower-/c waves, 
this mechanism allows waves generated by low-energy particles to be converted to 
waves resonant with high energies, and this will enhance the rate of stochastic accel- 
eration of the electrons, for example. Also the theory of diffusive shock acceleration 
(DSA) greatly benefits from an inverse cascade: it is the wave intensities at the 
lowest frequencies that determine the maximum energies produced in a shock in a 
given time, but as the upstream spectrum of plasma waves is usually generated by 
streaming instabilities due to the accelerated particles themselves, the spectrum of 
waves is usually an increasing function of frequency in the range resonant with the 
highest energy particles. 10 Obviously, an inverse cascade will increase the intensity 
of waves resonant with the highest energy particles and, therefore, increase the rate 
of acceleration. 

The previously analysed high-beta interactions, 

1+ ^ R+ + LT 
1+ ^ L+ + R- 
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are also extended to higher frequencies. Both of these reactions are also possible 
in low-beta plasmas, as long as the frequency of the L wave is close to the ion- 
cyclotron frequency. As sound waves are rapidly damped in general, the reactions 
usually proceed to the direction of wave coalescence. Thus, these reactions may 
represent a new dissipation channel to the R waves in low-beta plasmas, although 
cyclotron damping of the L waves is also large near the resonance. 

We have also discovered completely new interactions, which are not possible in 
an MHD description. In these reactions, all three waves propagate in the same 
direction, 

I + <-> R+ + L+ 
L+ <-> L+ + I+ (low (3) 
R + <-> R+ + I+ (high/3), 

which is strictly forbidden for non-dispersive waves. The latter two interactions 
shift energy from the dispersive frequency range to the MHD range, but an inverse 
cascade does not really develop, as the wave falls out of the frequency range able 
to decay already after the first interaction. This process, nevertheless, may help us 
understand, why spectrum of magnetic fluctuation in the solar wind experiences a 
clear break above the ion-cyclotron frequency, 11 although the right-handed mode 
should not experience cyclotron damping in this frequency range. 

Assuming that the sound waves are rapidly damped, all the discussed three- wave 
interactions provide a means of plasma heating. As the dissipation mechanism is 
not based on cyclotron resonance, the resulting heating rates on ions and electrons 
may be significantly different from direct ion-cyclotron damping. 

In this paper, we only made use of the symmetries of the response tensor to find 
out the interactions with non-zero interaction rates. In a follow-up paper, we will 
study numerically the interaction rates of the three-wave interactions presented in 
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this paper. This will allow us to deduce the importance of the interactions relative 
to other plasma phenomena, and determine the effects of three-wave interactions on 
plasma heating rates and particle acceleration. 

In future work, it will also be important to extend the study also to waves 
propagating at a finite angle with respect to the magnetic field. In this case, as the 
polarisation tensor will become more complicated and the response tensor loses its 
transverse isotropy, many reactions should be allowed to occur. Another potentially 
important point to include to the model is multiple ion species. As the L waves will 
develop more resonances and subsequent cutoffs, especially the interactions L — > 
I + + L ± will be different from the ones presented here. 

In conclusion, three-wave interactions in collisionless plasmas involving dispersive 
waves show a richer variety of possible reactions than the MHD counterparts. These 
reactions are potentially important in many branches of plasma astrophysics ranging 
from particle acceleration to plasma heating. We have identified some of them, but 
the theory has many more application not discussed in this paper. 
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